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We consider the spontaneous creation of a dc voltage across a strongly coupled semiconductor 
superlattice subjected to THz radiation. We show that the dc voltage may be approximately pro- 

■ portional either to an integer or to a half-integer multiple of the frequency of the applied ac field, 
(~| ' depending on the ratio of the characteristic scattering rates of conducting electrons. For the case of 

, an ac field frequency less than the characteristic scattering rates, we demonstrate the generation of 

^ ' an unquantized dc voltage. 
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The theoretical analysis of nonlinear transport properties of strongly-coupled semiconductor superlattices (SSLs) 
irradiated by a high-frequency electric field began already in the mid 1970's [1]. Recently, experimental progress 
in creating powerful sources of THz radiation, the development of a coupling technique [2,3], and improvement in 
I the fabrication technology of microstructures leading to very high carrier mobility [4] have stimulated many new 
O ' theoretical investigations of this long-studied problem. Among these works have been studies of strongly nonlinear 
O , , effects including multistability [5] , short pulse generation [6] , chaos [7-9] , and spontaneous generation of a dc voltage 
in a purely ac-driven SSL [10-13]. 
. In this paper we investigate the last effect in further detail and discuss the appearance of new dc voltage states, 
' which are a generalization of the integer dc voltage states in SSLs previously described [10,11]. These states are 
related to the complex dynamics of miniband electrons in an SSL and the formation of the Wannier-Stark ladder in 
a purely ac driven SSL. Two distinct mechanisms are known for the spontaneous generation of a dc bias in purely 
I ac-driven SSLs [10,11,13]. Both these nonlinear mechanisms work in SSLs with a high mobility and a relatively high 
' level of doping, when the effects of a self-consistent electric field generated by an electron's motion become significant. 
The first mechanism arises if the ac field frequency, cu, is much greater than the plasma frequency, ujpi, and is related 
[10,13] to an instability caused by absolute negative conductivity in the ac-driven SSL [14]. The effect has been 
attributed [10] to the phenomena of dynamical localization of electrons [15] and miniband collapse in a coUisionless 
SSL [16]. It was reported that the generated dc bias is such that the "induced Bloch frequency" ojb = eaEdc/fi [Edc is 
the spontaneously generated dc electric field and a is the SSL period) is approximately equal to the ac field frequency 
: ^ [10,13,17]. 

Ch . The second mechanism is responsible for dc bias generation when the ac field frequency is near the plasma resonance, 
^ ■ a; ~ Upi] it arises for a smaller ac field strength than is required in the previous situation [11]. The instability 
responsible for the dc bias generation in SSLs with strong enough electron scattering also results in chaotic motion 
in the case of small scattering rates or for a coUisionless SSL [7,11]. The creation of a dc bias may be qualitatively 
explained [12] and classified [13] using the semiclassical theory of wave- mixing in SSLs [18]. 
^ ' In this paper we re-examine the problem of spontaneous dc voltage generation in an SSL subjected to a THz 
5^ \ electric field. We show that, depending on the relative values of the scattering rates and the ac field frequency, a 
variety of different dc voltage states can exist, including both integer and half-integer quantized states, for which the 
induced Bloch frequency is approximately an integer or half-integer multiple of the ac field frequency, and completely 
unquantized states. In particular, if the electron velocity relaxation rate, 7„, is sufficiently different from the electron 
energy relaxation rate, 7^, and lo > 71,, we find integer states with lob ~ nuj (n = ±1, ±2, . . .); while for 7^ — 7^, 
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the states are close to the half-integer states ljb ~ nuj/2. In contrast, in the case of low-frequency driving or high 
damping, w < (7i;,7e), the dc voltage states are unquantized. 

We study electron transport through a single miniband, spatially homogeneous SSL with period a and miniband 
width A, which is subjected to an ac electric field E{t) = Eq coscut along the SSL axis. For the tight-binding cnergy- 
quasimomentum dispersion relation e{k) = (A/2) [1 — cos(A;a)] (A; is the electron wave vector along the axis of SSL), 
the dynamics of electrons is described by the superlattice balance equations [10,11,13] 

V = UW — 'JyV, 

= -UV - ^s{w - Weq), (1) 

il = WpiV - au + Iext{t), 

where v = moVa/h, w = {s — A/2)(A/2)~^ and Weq are a scaled electron velocity, a scaled electron energy, and an 
equilibrium value of scaled electron energy, respectively, and toq = (2?i^)/(Aa^) is the effective mass at the bottom of 
miniband. The scaled variables v{t) and w{t) are proportional to the variables V{t) and e{t), which are the electron 
velocity and energy averaged over the time-dependent distribution function satisfying the Boltzmann equation. The 
lower (upper) edge of the miniband corresponds to w = —1 (w = +1), and the value of Wpq is a function of the lattice 
temperature (for a thermal equilibrium Wgq < 0). The variable u{t) is related to the electric field inside the SSL E{t) 
as u = eaE/h. In deriving Eqs. (1), we assumed that the electrical properties of an SSL of total length I and cross- 
section S can be modeled by an equivalent high-quality circuit [10] which consists of a capacitor C = (eoS')/(47rZ) 
(eo is the average dielectric constant for the SSL) driven by an ac current of the form I^xt = — w^wsinwi, where 
u)s = eEoa/h, in our scaled units. 

The first Eq. of set (1) describes an acceleration of electrons under the action of electric field E{t) and their slowing 
down caused by an effective friction due to scattering. In original dimensional variables, the term uw in r.h.s. of the 
first Eq. looks like eE/m{e), where m{e) = mo/{l — 2e/A) is energy-dependent effective mass of the electrons in 
SSL's miniband. The second Eq. describes a balance of electron's energy gain under the action of electric field and 
energy loss due to scattering. Finally, the third equation describes a balance of difi'usive, external and displacement 
currents in the SSL. The degree of nonlinearity in Eqs. (1) is controlled by the value of miniband plasma frequency, 

1 /2 

Wpi = (47re^7V/moeo) , which is a function of the electron doping density N, while the parameter a determines the 
quality of the effective circuit (a ^ ujpi , oj) . 

The relaxation processes for miniband electrons are characterized by an average energy scattering rate, 7£, as well 
as by an average velocity scattering rate, 7^=76+ 'Jeh where 7ei is an average rate of elastic collisions [3,10]. The 
scattering rates, and 7^, can have different values depending on the material, the doping density, the temperature, 
etc. In particular, for microstructures with modulation doping [19], the ionized impurities are spatially separated 
from electrons, which greatly reduces the elastic scattering rate 7ei so that 7^ w je- In contrast, for many vertical 
SSLs operating at room temperature, the scattering rate for electron velocity is about of order of magnitude greater 
than the characteristic scattering rate of electron energy, 7„/7e « 10 [3,10]. 

We solve the nonlinear balance Eqs. (1) numerically for the initial conditions v{0) = 0, w{0) = Weq = —1, with the 
circuit damping rate a/ujpi = 0.01, and for two typical sets of relaxation constants: (i) jv/je = 10, js/<-^pi = 0.01, 
and (ii) "fv/je = Ij lel^pi = O-l- After removing the transients, we calculate the time-average value (u), which gives 
the value of the Bloch frequency determined by the spontaneously generated dc bias E^c.'- = eaE^d^ = (u). Figs. 
1 and 2 present the results of computations of (w) for 201 value of Ug equally distributed in the range < < 2upi 
for each driving frequency. 

For the first set of relaxation rates and for u > 27^, plateaus of integer-quantized states are clearly observable: 
lub ~ niu, with n = ±1,2,3,4 in Fig.l. However, at low frequencies, uj < 2"/y, instead of steps there is a region of 
unquantized dc states. The dependence of the induced Bloch frequency (u) on the ac field frequency lu for the second 
set of damping parameters is presented in Fig. 2. Here, qualitative differences from Fig. 1 appear, including 1) the 
existence of half-integer states, specifically states with n w ±1/2, and 2) the nonzero width of the n = plateau. One 
can expect the width of a plateau to be equal to the scattering rate, O.lujpi in this case. As a result, wc have plotted 
the fines n ± O.lujpi/u) for n = 0, ±1, ±2, ±3, —4, —5 in Fig. 2. We found that, indeed, most points for a given plateau 
lie in the region demarcated by the two lines with same n. In contrast to Fig. 1, there are states with {u)/uj < 0.1, 
but (u) ^ 0. Also, some points fall very near the line {u) = ±0.5a;. As an example we refer to the solution of Eqs. 
(1) for u/wpi = 0.6 and ojg/^pi = 0.6 (other parameters are same as in Fig. 2), which in the phase space corresponds 
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to a symmetry-broken limit cycle-'^ with (u) = 0.289ujpi. 

If u)/'y is less or order of unity but lext/^pi is large enough and 7„ = 7^ = 7, then the existence of nonquantized 
and half-integer dc voltage states can be demonstrated analytically. In these limits, the motion on an attractor of 
dynamical system (1) is governed by the following pendulum equation [20] 

e + je + {-wegpli sine = -2hxt{t), (2) 

where we made the substitutions: v = {Wf,q/2) sva.6, w = (weg/2)(l + cos^) and 

u = ^/2 + 7tan(6l/2). (3) 

The dc voltage (u) can be obtained as a result of averaging over time of r.h.s. of Eq. (3). The equation (2) looks same 
as a motion equation of the well-known Stewart-McCumbcr model from the theory of ac-drivon Joscphson juctions 
[21]. However, the principal difference from this model also exists. In the Stewart-McCumber model the voltage across 
the Josephson jmiction is proportional to the velocity of pendulum [21], while in our case the voltage is a function of 
both velocity and co-ordinate (see Eq. (3)). This difference plays a principal role in the explanation of existence in 
our case of both unquantized and only approximately quantized dc voltage states, which are absent in the Josephson 
junction model. 

The damped and driven pendulum (Eq. (2)) have two distinct types of attractors with regular dynamics: rotating 
and oscillating [21,22]. Majority of rotating states are phase-locked, i.e., {6) = {n/l)u) [n and / are integer numbers); 
they are observable mainly at w ^ 7 [22] . As it is evident from Eq. (3) , such phase- locked rotational pendulum states 
are responsible for the generation of half-integer dc voltage states in SSL: {u) ~ {n/2l)oj. Note that these corrections 
to these quantized states arise from the contribution of ^-dependent term in Eq. (3). The corrections can give the 
dependence of (u) on ac current amplitude, Wg, its frequency, w, and scattering constant, 7. However, it can be shown 
that a relative contribution into these corrections is controlled by the parameter j/iVpi [20] . Therefore, the dependence 
of approximately quantized voltage on ac field amplitude and frequency is weak until j/ujpi <^ 1. 

For oscillating attractors, the equality {0) = is always valid [22]. If additionally stationary pendulum oscillations 
are symmetric, i.e. {6) = 0, then the dc voltage can not be generated, {u) = 0. However, the pendulum (Eq. (2)) 
can demonstrate symmtry-broken oscillations, for which (6) ^ [22,23]. The symmetry-broken oscillations mainly 
exist at the low frequencies of external force; they can survive at w < 7 even for a strong damping, when j/tVpi ^ 1 
[24]. As it is evident from Eq. (3), just the symmetry-broken oscillations corresponding to 0) = and {0) ^ 0, are 
responsible for the generation of the unquantized dc voltage in SSL. 

Unquntized dc voltage states exist in the low frequency region for different rations of the scattering constants. In 
order to understand the transition from quantized to unquantized dc bias states as the scattering rates increase while 
maintaining 7„ 7^, we present in Fig. 3 the dependence of (u) on lu for strong damping. As it is evident from this 
figure, the strong damping destroys all quantized states; dc bias generation persists only for some unquantized states. 

It is instructive to consider the mechanism of unquantized dc bias generation in the terms of energy levels structure. 
We offer the following qualitative explanation of the underlying physics: when the electron scattering rates are 
sufficiently small and the amplitude of the ac field is large enough, the SSL spontaneously creates a Wannier-Stark 
ladder with the spacing, wb, that makes multiphoton absorption of ac field most effective, i.e. u>b ~ nw. In the 
pendulum analogy this means the appearance of the pendulum rotations, which frequency of rotation is proportional 
to the generated dc bias. However, for larger damping, the rotations are ceased and therefore the symmetry broken 
oscillations remain, that means that only a small bias can be generated in the SSL. Hence the spacing of the Wannier- 
Stark ladder states is less than the ac field frequency except for very low frequencies 7. In this case, the broadening 
of the self-organized ladder levels is quite comparable with their spacing; it is therefore practically impossible to 
achieve quantized values of the voltage (or a rotation of the pendulum). In a real space picture such a situation 
should correspond to an appcarcncc of some kind of semi-localized Wannier-Stark wave functions; this structure of 
the ladder is reminiscent of the wave function of biased short superlattice [25,26], where the spacing between energy 
levels depends non-linearly on the bias voltage [26]. In our case, a dc bias is created by ac field, thus we have a weak 
dependence of ladder spacing on ac field strength via the self-induced bias. 

We have performed a systematic numerical study of the positions and widths of different plateaus and of the 
unquantized states for many values of the driving amplitude and frequency and several different initial conditions at 



^This is a limit cycle whose projection on the v — u plane is not symmetric about the origin, in contrast to a symmetric limit 
cycle (see [11]). 
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different damping levels [27]: the results are in a qualitative agreement with situation described above. We now make 
several remarks on the results of this search. First, we found no indication of the noninteger (fractional) dc states 
for 7^ 7^ 7e- However, for 7^ = 7e) we find that the 1/2-dc-states are quite common. Moreover, for weak enough 
damping, we additionally saw a few dc states that are very close to fractional states of the form n/k with n being 
an integer and k always being an even integer. Such dc states are formed by the symmetry-broken limit cycles with 
large even periods. As examples, we refer to 3/2-states formed by period-12 and period- 24 limit cycles, as well as to 
the 7/6-statc formed by a period-12 limit cycle; both occur for 7t, = 7^ = 0.05cjp; and a./ijjpi = 0.01. We should note, 
however, that for weak damping, chaotic behaviour is quite typical [7,11] and that both the stable limit cycles and 
the asymmetric chaotic attractors, which are responsible for the generation of stable quantized dc voltage states in the 
SSL, occupy only a small amount of parameter space of the system [27,20]. 

Importantly, it appears possible to achieve the parameter values used in our simulations in the conditions of 
experiments. To begin with we refer to the experiments with a heavily doped SSL (N = 8 x 10^^ cm~^) having at 
room temperatures scattering constants jy w 10^"^ and 7^ ^ 0.l7t, [3]. Using the values a = 4.8 nm, A w 50 meV 
and eo ~ 13 [3], we get uipi = 1.2 x 10^^ rad/s and j/copi « 0.84 (compare, e.g., with data of Fig. 3). 

Longer scattering time provide SSLs based on the cleaved edge overgrown technique: 7"^ « 3 ps at low temperatures 
[4]. In conditions of experiment [4], for a = 10 nm, A « 20 meV, N = Ngb w 3 x 10^''' cm^'^ (Ng = 3 x 10^^ cm^^ is 
an electron areal density and b w lO"** cm is sample's thickness), the miniband plasma frequency is ujpi = 3.1 x 10^^ 
rad/s providing j/ojpi « 0.1 (c/. with data of Fig. 2). In our numerical simulations we used the ac field strengths 
satisfying lUs < 2tjp;, what in the physical units correspond to the realistic values Eq < 5 kV/cm. Necessary for the 
observation of our effects ac field frequencies belong to the THz range. Finally, we should also note that for all listed 
cases, the standard condition of validity of the single-miniband approximation, A ^ TiWg [14], is well satisfied because 

wsn/i\ < 0.1 . 

In summary, we have shown that a semiconductor superlattice irradiated by a high-frequency electric field can 
spontaneously generate a dc bias, which can be quantized in approximately integer or particular fractional ratios of 
the driving frequency, or completely unquantized. In this respect, the effects in semiconductor superlattices are no 
less rich than their counterparts in Josephson junctions subjected to a microwave field, where the exactly integer and 
the exactly fractional dc voltage states ("phase-locked states") are known [21]. 
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by the Academy of Finland (grant 163358) and NorFA. 
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FIG. 1. The dependence of spontaneously generated dc bias {u)/uj on ac frequency uj, scaled to the miniband plasma 
frequency ujpi, and for 7^ = O.lujpi, 7e = a = O.Olujpi. 




FIG. 2. Same as in Fig.l, but for 7^ = 75 = O.lcjpj, a = O.Oluipi. 



5 



0.1 




0.4 0.8 1.2 1.6 

CO 

FIG. 3. Same as in Figs.l and 2, but for strong damping: 7„ = u!pi, fe = O.lWj,;, a = 0.03Wp(. 
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